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Abstract 

Within Dyson-Maleev (DM) transformation and self-consistent mean-field 
treatment, the Neel order /disorder transition is studied for an antiferromag- 
netic Heisenberg model which is defined on a square lattice with a near- 
est neighbour exchange J\ and a next-nearest neighbour exchange J2 along 
only one of the diagonals. It is found that the Neel order may exist up to 
J2/J1 = 0.572, beyond its classically stable regime. This result qualitatively 
improves that from linear spin- wave theory based on Holstein-Primakoff trans- 
formation. 
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The two-dimensional (2D) antiferromagnetic Heisenberg models have attracted great in- 
terest in recent years, partly because of the fact that the parent compounds of the high 
temperature superconducting materials are excellent realizations of quasi-2D quantum anti- 
ferromagnets ji] . While the unfrustrated Heisenberg model has been well understood, much 
attention has been paid to the frustrated models such as square-lattice nearest and next- 
nearest neighbour interaction (so called J1-J2) model, triangular lattice model and Kagome 
lattice model etc In these frustrated Heisenberg models the property of the ground 

state, whether magnetically ordered or disordered, is a subject of considerable interest. For 
example, the J1-J2 model takes on Neel order at small J2/J1 and collinear order at large 
J2/J1, which are seperated by a region of disordered state 0. 

Very recently a S = 1/2 Heisenberg model, which defined on a square lattice with a 
nearest neighbour antiferromagnetic exchange Ji and a next-nearest neighbour exchange J2 
along only one of the diagonals of the lattice as shown in Fig. 1, has been proposed [§,0. 
Its Hamiltonian is written as 

H= J 1 ^S i -S i + J 2 £S,-S m , (1) 

{ij) (Im) 

where the notation (ij) denote nearest neighbour bonds and (Im) denote next-nearest neigh- 
bour bonds along only one diagonal. Topologically this model is equivalent to the Heisenberg 



ant if erromagnet on an anisotropic triangular lattice ||1 0| , |1 1|| . In special cases J2 — 0, J\ — J2 
and Ji = 0, it will recover to unfrustrated square lattice model, isotropic triangular lattice 
model and decoupled spin chains, respectively. Therefore this model provides a way of inter- 
polating between several well-known one and two-dimensional models. One can study the 
role of frustration in going from one to two dimensions. On the other hand, this model is of 
direct relevance to the magnetic phases of some quasi-2D organic superconductors. It has 
been argued that this model may describe the spin degree of freedom of the insulating phase 
of the layered molecular crystals k-(BEDT-TTF)2X |12|. The parameter J2/J1 for these 
materials is suggested to be ~ 0.3 — 1 and the magnetic frustration will play an important 
role. 
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FIG. 1. The 2D square lattice with nearest neighbour interaction ,J\ and next-nearest neighbour 
inteaction J2 along only one of the diagonals. 

Classically, the ground state of the model can be derived straightforwardly as a function 
of the ratio J%j J\ if we assume that the spins lie in the xz plane and are described by a spiral 
form Si = S^sin^j, 0, cos^j) as shown in Fig. [l|. Here the angle 9i = qr, and the wavevector 
Q = (OS q) defines a relative orientation of the spins. Minimization of the classical energy 
with respect to q gives the result that the ground state take on Neel order (i.e., q = ir) for 
J2/J1 < 1/2, and spiral order with q = arccos(—Ji/2J 2 ) for J2/J1 > 1/2. 

The quantum model has also been studied numerically and analytically. The series 
expansions were adopted by Zheng et al. [Q for numerical calculation. It was found that the 
Neel order persists up to J2/J1 = 0.7. In the region 0.7 < J2/J1 < 0.9 there is no magnetic 
order and for larger values of Jij J\ there is incommensurate or spiral order. It is interesting 
to note that the Neel order exists beyond its classical result {J%j J\ = 0.5). Analytically 
the standard and simple linear spin-wave theory (LSW) based on Holstein-Primakoff (HP) 
transformation was used by Merino et al. to discuss the possible ordered and disordered 
states. It is helpful to repeat some details here for later discussion. First, for convenience 
the spin at each site is rotated along its reference direction characterized by the angle 9i. 
Such a rotation may be accomplished by the following transformation for spin operators 

S? = sin0 i 5? + cos0 < S? 
cv _ qy 

J i — J i 

S* = cos 6$ -sin 6iSf 

For new spins S« the reference ground state becomes ferromagnetic. The rotated Hamiltonian 
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becomes 

H = Ji YfcoaOiASfSj + SfS?) + sin % {S*S* - S*S*) + 5fSj] 

(ij) 

+j 2 J2 [cos M^f <fe + s?s* m ) + sin MSfs* - s?s*j + sfsyj . (2) 

(lm) 

with 6^ = 6i — Oj = q, Q\ m = 9\ — 9 rn = 2q. Then with HP transformation the above 
Hamiltonian with only quadratic terms kept may be diagonalized in the momentum space. 
The final dispersion relation for the spin excitation is 

^ = ^{[J(k + q) + J(k - q)]/2 - J(q)}[J(k) - J(q)] 

with J(k) = Ji(cosk x + cosk y ) + J 2 cos(k x + k y ). By calculation of the magnetization (S?) 
(see also the dashed line in Fig. ^ later) it was suggested by the authors in Ref. that a 
possible Neel order /disorder transition happens at the point J 2 j J\ — 0.5. However, we want 
to point out the inherent limitation for LSW theory here. From the spectrum u^, it is easy 
to find that, in order to ensure the argument in the square root always positive in the whole 
Brillouin zone, the parameter q has to be set as ir (characterizing Neel order) in the region 
J2/J1 < 0.5, but must not be set as it in the region J2/J1 > 0.5. This means that within 
LSW theory the Neel ordered state can never appear beyond its classically stable region, 
which is in contrast with the numerical result || . As also pointed out by the authors in Ref. 
themselves, the interaction between spin waves becomes very large at the transition point 
and it may lead to a completely different picture for the states. So it is very necessary to go 
beyond the LSW theory to see how the above result will be modified, which is the purpose 
of this paper. Instead of HP transformation, the Dyson-Maleev (DM) transformation will 
be adopted to avoid 1/S expansion. It has been also recognized that the DM transformation 
must be prefered if one needs really to go beyond LSW theory within a perturbation scheme 
||13|| . In the treatment of the so called J\-J 2 model, it has been shown that the spin-wave 
theory based on DM transformation gives perfectly consistent result as that from numerical 
calculation; much better than that from LSW theory jj. 

In this work we will not discuss the possible spiral state for large J 2 / Ji, but focus 
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on the Neel order /disorder transition in the intermediate J2/J1 region, especially on the 
problem whether the Neel order may appear beyond its classical region of stability which is 
one of the most interesting topics for this model. Technically, the Hamiltonian under DM 
transformation have terms as high as sixth order when spiral state is considered (i.e., for 
general q), which are relatively complicated to treat. When only Neel state is considered the 
transformed Hamiltonian has no term higher than fourth order, which can be easily treated 
by mean-field (MF) theory or perturbation theory. Explicitly, one may apply the DM 
transformation onto the original Hamiltonian (|l|); or begin with the rotated Hamiltonian 
(§) by setting 9ij = it, 6i m = 2tt and then use the DM transformation for A and B two 
sublattices in the following form: 

= (1 — a\ai)ai, S~ = aj, S- = 1/2 — ajaj 
Sr = b]{l - b]bj), S+ = bj, S* = 1/2 - b]bj , 

where a(a^), b(b^) are bosonic operators for sublattices A and B, respectively. Then the 
Hamiltonian (0) is transformed into 

H = -J x £(1/2 - oJoi)(V2 - b)bj) + [(1 - aJoOai&j + a\b](l - b]bj)]/2 

(«> 

+ J A H (1/2 - ajai) (1/2 - a^arn) + [(1 - a]ai)aial n + (1 - a^a^arnaj}/ 2 + £ a -> b} . (3) 

(lm)eA (lm)eB 

To diagonalize the above Hamiltonian, we treat the quartic terms with a self-consistent MF 
theory. For those terms which could be decoupled in two ways, we will combine the two 
kinds of decoupling form together through a weight factor A as introduced by Chu and Shen 
nj. For example, the term like a\aib^bj will be decoupled in the way 



ajaib^bj ~ \[{a\ai)b^bj+a\ai{b^bj) —{a\a^ ^ . 
The parameter < A < 1 reflects the competition between two decoupling ways. Its value 



may be decided by minimization of the energy, which was found to be 1/2 in Ref. for 
unfrustrated lattice; or it is required to be equal to 1/2 in order to keep the symmetry of the 
Hamiltonian before and after decoupling ||15|| . We will take the value A = 1/2 throughout 



our calculations. With definition of several parameters: u = (a|oj) = (fef&i), u = (affej) = 
(ciibj), w = (a\a m ) = (b]b m ), we may obtain a quadratic Hamiltonian in the momentum 
space: 

H = Y}°^ a W + b %) + E k (a k b^ + 4&U)] + const > ( 4 ) 

k 

C k = 2Ji(l — u + v ) — J 2 (l — w + w)[l — cos{k x + fey)] , 
£7 k = — Ji(l — it + v)(cos fe x + cos fey) , 
const = — [J x (-1 + 2m 2 + 2v 2 - Auv) + J 2 (l - 2m 2 - 2w 2 + 4utw)/2] , 

where the summation is over half of the original Brillouin zone and iV is total number of 
lattice sites. Under a Bogoliubov transformation 

a k = cosh Ak ak + sinh Ak 6_ k 
&L k = sinh Ak a k + cosh Ak &l k 

with tanh2Ak = — -Ek/Ck, the Hamiltonian (f|) may be diagonalized into 

# = E^k(4«k + b[b* + 1) - C k ] + const (5) 



with the excitation spectrum cu k = yC k — -E^- Correspondingly the self-consistent equations 
for u, v, w are expressed as 

1 x — > Ck / 

u = — V . - 1/2 , 

1 ^ E k cosk x 
f = > 



,/C 2 ^ 2 



k ~~ -^k 



1 CkCOsffea; + fe„) 



The magnetization is simply given by m = 1/2 - u and the ground state energy is Eq 
Ek(^k - C k ) + const. 
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FIG. 2. The magnetization m calculated from our treatment (solid line) and that from LSW 
theory (dashed line) as a function of J2/J1. It goes to zero at about J2/J1 = 0.499 within LSW 
theory, but up to J2/J1 = 0.572 within our treatment. 




FIG. 3. The self-consistent results for the parameters v and w. 

We show the numerical results for the above self-consistent equations in Figs. |2] and |||. 
The magnetization m (derived from the parameter u) as a function of J2/J1 is given by the 
solid line in Fig. 2, which is the main result in this paper. It is found that the magnetization 
does not vanish until J%j J\ — 0.572, which gives the Neel order /disorder transition point. 
As comparison, the result from LSW theory is plotted by the dashed line; the magnetization 



goes to zero immediately before the classical value J2/ J\ = 0.5, see also Ref. [IT]. As 
expected, the current result is closer to the numerical one from series expansions and most 
importantly, it qualitatively improves the result from LSW theory. As we discussed before, 
the LSW theory is impossible to deduce a Neel ordered state beyond J2/J1 = 0.5. In our 
treatment the interaction between spin waves is actually partly considered, then the fact that 
Neel order may exist beyond its classically stable region is shown. It is quite possible that 
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the transition point will shift to larger value if the residual interaction between spin waves 
is included. This may be also hinted from Fig. ^ where the parameter v, which represents 
the antiferromagnetic correlation between the two original nearest-neighbour spins, is large 
near the transition point. 

For completeness, the ground state energy in the region < J 2 / ' J\ < 0.57 is also shown 
in Fig. § by the solid line, which is close to the result from LSW theory (the dashed line). In 
the small J2/J1 region the energy calculated here is a little lower than that of LSW theory, 
and becomes a little higher with increase of J2I J\- For J2/J1 > 0.5 the energy within LSW 
theory is derived from spiral state; there is a cusp at point J2/J1 = 0.5 |J. In the current 
case the state is still Neel ordered and the energy changes smoothly. 

-0.54 
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-0.62 
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J2/J1 

FIG. 4. The ground state energy in unit of NJ\ calculated from our treatment (solid line) and 
that from LSW theory (dashed line) as a function of J2/J1 (see text). 

In summary, within Dyson-Maleev (DM) transformation and self-consistent mean-field 
treatment, we have studied an antiferromagnetic Heisenberg model on a square lattice which 
includes a nearest neighbour exchange J\ and a next-nearest neighbour exchange J2 along 
only one of the diagonals. This model should be of direct relevance to some layered organic 
superconductors. In this work we focus on the discussion of Neel order /disorder transition 
for not large J2I J\- It is found that the Neel order may exist up to J2/J1 = 0.572, beyond 
its classically stable regime. This property is consistent with numerical finding from series 
expansions, which is one of the most interesting features for this model. Especially, because 
the interaction between spin waves is partly considered in our treatment, the result derived 
here qualitatively improves that from LSW theory based on Holstein-Primakoff transforma- 
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tion. It is certainly necessary to continue this work to study the large J2/J1 region where a 
spiral order will appear, so that a whole phase diagram may be constructed. 

The author would like to thank J. Merino for helpful discussion. This work was supported 
in part by Chinese NSF. 
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